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Abstract. We investigate the tunnelling zone V — m < E < V +m for a one-dimensional 
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potential within the Dirac equation. We find the appearance of superluminal transit times 
akin to the Hartman effect. 
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m: I. INTRODUCTION 

. In several recent articles, we have investigated in some detail one-dimensional electrostatic potentials 

by means of both the Schrddinger [1-3] and the Dirac equation [4-6]. Several original phenomena 
have been observed, such as the transition from resonance phenomena to multiple (infinite) peak 
formation [5] equivalent to a shift from wave-like to particle behavior in the barrier diffusion zone 
E > V a + 77i, where V is the barrier height, E one of the wave packet energies and m the particle 
mass. We have also investigated the compatibility of the barrier results with the Klein paradox [7-9], 
when m < E < V — m. In this latter case, we have noted the existence of dynamic localized states 
with a continuous spectrum [6]. These states are the nearest approximation to the bound states of 
■ Schrddinger or of Dirac in the evanescent energy zone considered in this paper. 

The evanescent zone is the last energy zone we need to consider to complete our analysis. It is 
given by V — m < E < V + m (V > 2m) or m < E < V + m (V < 2m). It has evanescent (real 
exponential) space forms in the barrier region. For a well potential, it is just such forms that give 
rise to discrete bound states. In this paper, we shall concentrate upon the single barrier potential 
and hence complete our analysis for this elementary structure. The evanescent stationary solutions 
become dynamic if instead of plane-waves we work with incoming wave packets. Then, the particles 
within the classically forbidden region are measurable only for a finite time, the time of transition from 
an incoming wave packet to reflected/transmitted wave packets. Even for the step potential in this 
energy zone there will exist, during this transitory time, current flow both into and then subsequently 
out of the step. Since, the stationary solution has zero current flow within the step, this feature is 
not always recognized. It is however obvious when one admits that there will be a non zero transitory 
dp/d/j for any space interval within the step. 

The most important barrier feature (both theoretically and experimentally) of this energy zone 
is tunnelling. A part of the incoming wave packet will continue its course beyond the barrier region. 
Its magnitude will be modulated by the barrier. For "large" barriers (compared to the wave packet 
size) an exponential reduction in amplitude cx exp[— ql] occurs, where I is the barrier length and q 
is y m 2 — (E — V,-,) 1 . This not only reduces the amplitude but modifies the transition spectrum. The 
smaller the transition amplitude is, the smaller the modifications in the reflected wave packet from the 
incident wave packet. However, in general, for both for the reflected and transmitted wave packets we 
have maxima in configuration space, and can apply the stationary phase method (SPM) to calculate 
reflection time delays and transition times [10]. With the Schrddinger equation the conclusion that the 
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transition time is independent of barrier width I, when I — » oo, is known as the Hartman effect [11]. 
Such a result is hard to avoid and, if the same occurs for the Dirac equation (subject matter of this 
paper and previously discussed by other authors [12-14]), we would have to face the unpalatable 
feature of superluminal velocities within the barrier. We warn that more than one type of transition 
time has been defined in the literature [10, 15-17] and for details we refer to reference [18, 19]. In this 
paper, wc intend to investigate this particular aspect of tunnelling by means of the SPM neglecting 
the possible ambiguities that this approximation is known to have. 

In the next section, we define all quantities and equations used. Some of these have been given 
also elsewhere [5, 6] but for completeness we re-present them. We also solve the stationary plane 
wave problem for the step and barrier. In Section III, we calculate the transition times by using the 
SPM and, based on numerical calculations, we discuss the appearance of superluminal velocities. Our 
conclusion are drawn in Section IV. 



II FORMALISM AND SOLUTIONS FOR THE BARRIER 

We shall work with a one-dimensional (electrostatic) potential in the Dirac equation. The chosen axis 
is the z-axis. However, we shall use the solutions and hence formalism of the full three dimensional 
case [21]. Thus, the stationary Dirac equation reads 

[E-V(z)} lo i){z)+i l3 ^{z)=mi){z) , (1) 

where 7 and 7 3 are two of the Dirac matrices (see below) and tp'(z) = dip(z)/dz. The explicit time 
dependence exp[— iEt] has been dropped from the above equation and hence ij) is only a function of 
the z coordinate. Our chosen representation for the gamma matrices is the Pauli-Dirac one, so that 



^={0 -l) ■ T=[-Vo) • (2) 

The barrier potential is fixed at V in the region < z < I and is zero elsewhere. The z axis is divided 
into three regions. Region I is the region of the incident and reflected waves (z < 0). Region II is the 
barrier region. Region III is that containing the transmitted wave (z > I). For V > m, it is convenient 
to divide the tunnel energy zone into two sub-zones both evanescent: (A) V < E < V + m and (B) 
V„ — m < E < V Q (V > 2m) or m < E < V (m < V < 2m). For V < m, only the evanescent zone 
(A), m < E < Vo+m, exists. The A zone corresponds to energy above potential (E > V ) and we will 
use the "positive energy" uW solutions modified for a non zero potential V a . The B zone corresponds 
to the below potential zone (E < V ). 

One of the questions we pose in this work is if all below potential solutions (" negative energy" ) 
represent physical antiparticles, be they oscillatory (free) or evanescent. For the B zone we shall use 
the it( s+2 ) solutions modified to allow for a constant non zero potential V . 

Spin flip is absent in all these problems (independent of the value of E) so by choosing an incoming 
spin up state, 

U M(p,E) = [l, 0, p/(E + m), Of, (3) 
wc find the following spinors in region II, 

A-zone: u^(±iq,E-V ) =[1, 0, ±iq/(E -V +m) , 0]*, ( ) 

B-zone: u^(±iq , \E - V \) = [ =R?/(|-E - V \ + m) , 0,1,0]*. U 

Only the spinor (p , E) appears in region III. 



A-ZONE: V < E < V + m (V > m) OR m < E < V + m (V < m) 

The solutions in the three regions are: 

Region I: z < 0, u^>(p,E) exp[ipz] + iJ> u^\— p , E) exp[— ipz] , 

Region II: 0<z<l, A > u^\iq , E - V a ) cxp[-gz] + B> u {l) (-iq , E - V ) exp[qz] , (5) 
Region III: l<z, T> w (l) (p,_E) cxp[ipz] , 
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and we are using un-normalized solutions but such that |i?>| is the reflection probability. The sym- 
bolism i?g (Tg) will be used for the A/B energy zones because E «s V respectively. The solutions 
in region II are the evanescent solutions oc exp[±gz]. We shall in what follows refer to the case of a 
step potential with only two regions (I and II), without treating this case separately we merely note 
that it corresponds to the above solutions with = 0. It should also be obtainable from the barrier 
solution when I — > oo, although some care must be taken when multiple peaks occur such as in the 
case of above potential diffusion [5]. The first of these barrier peaks reproduces the single step peak 
in the I — ► oo limit. 

Solving the continuity equations, ipi(0) = ipn(0) and ipu(l) — ipm{l), in matrix form yields 



1 1 

1 -1 



1 

R> 



where 



B> 



a 



and 



e -ql e ql \ 








ae~ ql -ae' 1 J 


. B > . 




T > . 



(6) 



E + m 



p E — V„ + to 

Since Dirac is a first order equation only continuity of ifi is imposed. Solving the above equations gives 



Thus, 

R> 



i 

R> 



l-a 2 
2a 



1 1 

1 -1 



1 1 

a —a 



-ql 



a ql 



a e 



-ae' 



r 









cosh(ql) + a sinh(gi) cosh(g/) 
cosh(ql) — a s'mh(ql) — cosh(ql) - 



smh(ql) T> cxp[ipl] and T> = cxp[— ipl) / 



sinh(g/) 
a 

sinh(g/) 



cosh(ql) 



,,<l>l 





' 1 " 




I— 1 



T e ipl 



2a 



sinh(gZ) 



(7) 



(8) 



The non relativistic limit, E — m = E NR <C m and V <JC to, reproduces the Schrodingcr results for 
the reflection and transmission coefficients (see the appendix for a detailed derivation). 



B-ZONE: V - m < E < V (V > 2m) OR m < E < V (m < V < 2m) 

For this zone, in the potential region, we shall use the spinor . Thus, we have 

Region I: z < 0, u^'(j),E) exp[ipz] + R < vS 1 >{—p,E) exp[— ipz] , 

Region II: 0<z<l, A< u^(iq , \E - V \) expf-gz] + B< u^(-iq , \E - V \) cxp[gz] , (9) 
Region III: l<z, T< u^\p , E) exp[ipz] , 



Continuity equations in matrix form yields 



1 1 

1 -1 



and 



where 



E+rn 
P 



1 

R< 

p-qi P 



E+m 





-P(3 




■(■ 


1 1 


) 








B< 







A< 
B< 



<1P 



(E + m)(\E-V \+m) 
Solving the above matrix equations, we find 



(10) 
(11) 



1 

R< 



l l 

l -l 



-p (3 
1 1 



-(3e~ ql f3e 
p-qi 



ql 

a ql 



cosh(g/) - 2i£M2i) 

3 



V 









cosh(gZ) 



cosh(ql) — (3 sinh(g/) 
— cosh(ql) — (3 sinh(g/) 



Jpl 





' 1 " 




1 



T >e l 



(12) 
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Thus, 



R* = 



1-/3 2 



sinh(qZ) T < exp[ipl] and T< =cxp[— ipl] / 



cosh(gi) — 



1 



20 



sinh(gZ) 



(13) 



Although not obvious, simple algebraic calculations show that i?< and T < are functionally identical 
to i?> and T> respectively (although, of course, valid in disjoint energy zones). Hence in the following, 
we will drop the suffixes and use 



T = cxp I — ip I + i arctan 
mV 



E -to 2 



EV 



'IP 



tanh(ql) 



i 



cosh 2 (ql) 



h — to 



EV a 



'IP 



sinh(ql) 



R = -i 



gp 



smh(ql) T exp[ipl] 



(14) 



III. TIME ANALYSIS 

In this Section, we shall calculate the analytic expression of the transition times by using the SPM 
and, by numerical calculation, we present the appearance of superluminal velocities. 



• SPM TRANSITION TIMES 



For this analysis the essential ingredient is the phase of the transmitted amplitude T . The gaussian 
envelope function g(p) will be assumed real and the wave packet function ^(x,i), defined in the 
standard way, is given by 

(15) 



\E»(a;,t) = J &pg{p)ip(z) exp[— iEt)] . 

A common choice for g(p) (unnormalized) peaked at p is 

g(p) = exp[-a 2 (p - p )7 4 ] , 

so that for the incoming wave Tp IAnc (z) = u^^p^E) cxp[ipz] the wave packet width is just a (large 
barriers are thus defined by 1/ a 3> 1). Due to the real nature of g(p), the phase in ^ , J)inc (a;, t) is simply 



0/.inc = P Z — E t . 

The SPM then sets the maximum of the incident wave packet at time t at 



(16) 



AE 



dp 



t 



Pa_ 

E 



t 



where E is the energy corresponding to the peak momentum value of p . Whence the incoming wave 
packet maximum reaches (ignoring interference effects with the reflected wave) z = at time t — 0. 

The SPM calculation of the transmission time uses the phase factor of T obtained in the previous 
Section, 

— tanh(gZ) 



arctan 



qp 



+ p{z-l) - Et 



(17) 



Deriving <f> UIjtTa . with respect to E and setting z = I, we find the following functional 



2E-V a 



qp 



E - m 2 - EV 



IP 



tanh(gi) 



{E -to 2 - EV )(V - E)l 
q 2 p cosh 2 (ql) 



lb 



E — to 2 — EV a 
qp 



tanh(gZ) 



(18) 
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The exit time of a single transmitted wave packet is then given by t tIB ,(E , 1), where E = \Jp 2 Q + m 2 
and p a is the peak momentum of the transmitted wave packet, i.e. (neglecting the spinor momentum 
dependencies) the maximum of g(p) \T\ (see Fig.l and Fig. 3). If the barrier is short then p w p 
because the gaussian dominates over the transmission amplitude. On the other hand, if we take 

I — » oo the functional t tIB ,(E, I) greatly simplifies giving 

tUE) ~ km tUE, = ^— - , (19) 
i^oo qp 

which is independent of I. At first glance, this would mean unlimited tunnelling velocities. Actually, 
this is often argued without taking into account the difference between p and p . The exit time is 
given by t(E ) and not by t(E ). If for example we set E to its maximum allowed value (compatible 
with tunnelling) E a = V a + m, we find 

T t rJ\V + m] — > oo , 

so that the tunnelling velocities are not in general unlimited. However, since we shall find in the next 
Section superluminal velocities for finite I, we shall not dwell upon the asymptotic tunnelling velocity. 
We only note that if instead of a gaussian wave packet (which technically overshoots the tunnelling 
zone) we use a truncated gaussian, we can avoid infinite SPM tunnelling times by truncating below 
V„+m. We warn however that a truncation in the momentum spectrum of a wave packet automatically 
introduces infinite wave packets in configuration space and we have to take care in using the SPM [2,5]. 

• SUPERLUMINAL VELOCITIES 

In the previous sub-Section, we derived an expression for the transmission time, 

Ur*(E a ,l) . 

This expression requires the knowledge of E , the peak value of the transmitted momentum distribu- 
tion. This implicitly assumes a single maximum. So, the SPM is certainly valid for moderate values 
of /, where the transmitted spectrum is almost gaussian as it is shown in Fig.l and Fig. 3 (I < 2a). 
In such a context, we have numerically calculated E (l) and whence obtained t tIB .(E ,l) by Eq. fTg)) . 
In Fig.2-a and Fig.4-a, we have plotted t tla (E ,l)/a against I /a. We observe that for i > a, we have 
to use instead of the maximum momentum value, the average value of the spectrum, i.e. E a = (E). 
This is because for a momentum curve which ends upon a maximum, and it is thus very asymmetric, 
model numerical calculations show that a more accurate result for SPM times is obtained with the 
use of (E). The surprising feature of the curves given in Fig.2-a and Fig.4-a is obtained by taking the 
ratio of the coordinates, 

t t r a (E ,l) 

the effective velocity of the tunnelling processes. This is plotted in Fig.2-b and Fig.4-b as function of 

I I a. As it can be seen from Fig.2-a and Fig.4-a, there is a plateau region where the transit time is 
independent of the value of I /a. In this region the effective velocity grows linearly (see Fig.2-b and 
Fig.4-b). However, the surprising feature is the numerical value of the velocity in this region: it is 
already greater than one both for the relativistic and for the non relativistic case. We have no need 
to go to the infinite barrier width limit (Hartman effect) to find superluminal velocities. 

IV. CONCLUSIONS 

We have studied in this paper the tunnelling phenomena predicted by the Dirac equation. One of 
the principle questions posed at the start was if a Hartman effect exists also for the Dirac equation. 
The answer is positive since the spinors play no significant role in the calculation of the transmission 
times. There is a difficulty with the fact that the different momentum dependence of the spinors 
lead to different transmission times for the components. However, this does not modify the result of 
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each exhibiting a Hartman-like effect. It is the SPM which obliges us to work with wave function 
amplitudes rather than with the transmission probability function, in which spinor components have 
been summed over. However, if all spinor components yield supcrluminal velocities then superluminal 
velocities must be expected. The Hartman limit (I 3> a) has an added complication because the 
transmission function dominates the incoming convolution function (gaussian in all our calculations) 
and it is not yet clear how, or even if, the SPM works in the absence of a clean maximum in the 
momentum distribution. We have avoided entering into this equation because there is no need to go 
to I S> a in order to exceed the velocity of light. 

A second question involved in our study was the identification of the nature (charge) of the 
particles temporarily (for wave packets) in the classically forbidden barrier region. This is a relevant 
question when one recalls that in the Klein energy zone (E < V — m) antiparticles are created 
and/or annihilated in the barrier region. For an antiparticle the barrier becomes a well and the 
mathematical below potential (V ) particles of energy E are re-interpreted as physical above potential 
(— V ) antiparticles of energy — E [6,9]. It is tempting to consider all "particles" with E <V (below 
potential) to be in fact physical antiparticles, even if associated to evanescent terms (B tunnelling 
zone). This is the reason, we divided the tunnelling energy zone into two the A ( Vq < E < V„ + m) 
and the B (V — m < E < V ) zones. We now give an argument based upon our studies of the Klein 
zone that says that this hypothesis is not true. 

Let us consider in the following the simple step potential. The Dirac equation conserves proba- 
bilities. How is this consistent with pair creation in which more particles are reflected R > 1 (Klein 
paradox) than are incident? Physically, total charge is conserved but probability certainly is not. The 
answer is suggested by the well know fact that the below potential particles in the Klein zone have the 
"wrong" group velocity. This fact incidentally concords with the Feynman-Stuckelber conclusion that 
such below potential particles must actually travel backwards in time. Returning to our conundrum, 
we observe that in any formal numerical calculation, which ignores the antiparticle re-interpretation, 
the Klein paradox appears mathematically as the particular solution to a problem in which at t = — oo 
we have two opposite moving wave packets: the incident wave packet at z = — oo and the below po- 
tential wave packet (of appropriate size) at z = +00. When these two meet at time t ~ at the step 
discontinuity z = the continuity equations tell us that they unite and form a single wave packet, 
the reflected wave packet [20]. In this way probability is indeed seen to be conserved. It is the re- 
interpretation in physical particle/ antiparticle terms which alters our viewpoint. However, the Dirac 
equation (with a real potential) can be viewed in this mathematical picture with only particles, albeit 
with energies both above and below potential. Indeed this is the way everyone treats the stationary 
plane wave problem, including ourselves [6] . 

Let us now apply the same mathematical viewpoint to the step in the energy zones A and B. In 
this case there cannot be any effective particle flow from z = 00 since the stationary solution under 
the barrier is, in both energy zones, a pure exponential decreasing space function. At time t = —00 
only the incoming wave packet exists in this case. Eventually, for t = +00 only the reflected wave 
packet exists. For times within the transmission period during which complete reflection occurs, we 
cannot have a reflection coefficient R > 1, even if only for an instant, without violating probability 
conservation. This conclusion is independent of the choice of the A or B zones. It is based upon the 
impossibility of having a modification of the initial conditions so as to reproduce a transitory Klein- like 
paradox. Consequently, even in tunnelling, the probability density under the potential must represent 
the same particles as those of the incoming wave. 

It would be desirable to conclude the debate on superluminal velocities in tunnelling phenomena 
but, at the moment, the results are far from being conclusive. For a potential of the order of the mass, 
the Dirac equation cannot be viewed as a one-particle equation, and particle creation is expected to 
play an important role. There is still much to be studied in potential problems within field theory and 
the final answer can only be reached by analyzing tunnelling phenomena in a second quantized the- 
ory [21,22]. However, this topic exceeds the scope of this paper and it will be appropriately discussed 
in a forthcoming article. In such a spirit, this paper has to be seen as an initial work to stimulate 
further investigations. 
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APPENDIX 



Let us discuss in detail the non relativistic (NR) limit, E — m — E NR <C m and V a <C m. We recall 
that for V < m only the evanescent zone (A), < E NR < V , exists. For the convenience of the reader, 
we rewrite the Dirac reflection and transmission coefficients given in the text, 



T = cxp[ — ipl] I 
1-a 2 



cosh(ql) 



1 + a 2 



2a 



sinh(qZ) 



R 



2a 



smh(ql) exp[ipZ] , 



where 



p = \/ E 2 — m 2 , 

q = ^m 2 -{E~V f 

q E + m 
a = i — 



p E — V + rn 



Taking the NR limit, wc obtain 



V ->• Pnr = \J 2m E KR , 

q -> <?nr = \/2m(V - £ m ) 



a — > i q NR I p NR 



Consequently, 



T NR = cxp[-ip NR l] j 
V 



cosh(q NR I) — i 



2-Enr — V a 



ii NR — i ■ 



2 v -E , N r(K — -Enr) 
The square modulus of these coefficients, 

2"nr 



2 y -Enr^K — E NR ) 
sinh((? NR T N R exp[jp NR Z] . 

4:E NR (V — E NR ) 



sinh(g NR I) 



|-Rn R | 2 — 



4 E NR (V - E NR ) + V 2 sinh 2 \yj2m (V - E NR ) I 

v: 



sinh" 



y/2m (V„ - E NR ) I 



(20) 



(21) 



(22) 



(23) 



4:E NR (V — E NR ) 

is often encountered in standard quantum mechanics textbooks (see for example [23] pag. 73). 



(24) 
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g(p) \t\ 
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Fig. i 
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13 



15 



ap 



Fig. 1. The transmitted momentum distribution is plotted as function of ap for different values of l/a, where 
a is the width of the incident wave packet at t = 0. The potential is equal to the mass of the particle, 
aVo = am = 10, and the peak of the incident momentum distribution is chosen to coincide with the center 
of the allowed zone (compatible with tunnelling) for the momentum, ap = a^V (Vo + 2m)/2 = 5\/3- For 
moderate values of l/a the transmitted momentum distribution is almost gaussian. The amplifications show 
the attenuation (due to the evanescent waves) of the transmission probability for increasing values of l/a. 
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Fig. 2-a 




o 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 l/a 



Fig. 2. This figure contains two curves. The plot in Fig. 2-a represents the variation of the (adimensional) 
transmission time, t trei (E , I) /a as function of l/a. The plot in Fig.2-b, represents the ration between the barrier 
width, I, and the transmission time Uib,(Eo, I), i.e. the effective velocity of the tunnelling process. The potential 
is equal to the mass of the particle, aV = am — 10, and the peak of the incident momentum distribution 
is ap — 5\/3- The maximum value of l/a has been chosen to be 2.0 in order to have an almost gaussian 
transmitted wave packet. This guarantees the validity of the SPM. The surprising feature of our numerical 
analysis is that the tunnelling velocity is already greater than one for a moderate barrier width. 
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Fig. 3. The transmitted momentum distribution is plotted as function of ap for different values of l/a. 
The potential aV a — 10 is now much smaller than the mass of the particle, am — 10 3 , and the peak of 
the incident momentum distribution is yet chosen to coincide with the center of the allowed zone for the 
momentum. ap = am^/Vo(V + 2m)/2 = 5^201 This case represents the non relativistic limit, V <m and 
E — m = _B NR <g m. The amplifications show a considerable attenuation of the transmission probability with 
respect to the relativistic case. 
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Fig. 4. The plot in Fig. 4-a represents the variation of the (adimensional) transmission time, t tIB .(E ,l) /a as 
function of l/a. The plot in Fig.4-b, represents the ration between the barrier width, I, and the transmission 
time ttm{Eo,l), i.e. the effective velocity of the tunnelling process. The potential aVo = 10 is much smaller 
than the mass of the particle, am — 10 3 , and the peak of the incident momentum distribution is ap = 5\/201. 
The maximum value of l/a has been chosen to be 2.0 in order to have an almost gaussian transmitted wave 
packet. This guarantees the validity of the SPM. The surprising feature of our numerical analysis, i.e. the 
tunnelling velocity greater than one for a moderate barrier width, is confirmed in the non relativistic limit. 



